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Abstract 

The solid-on-solid (SOS) model of an interface separating two phases is exactly sol- 
uble in two dimensions (d=2) when the interface becomes a one-dimensional string. The 
exact solution in terms of the transfer matrix is recalled and the density-density corre- 
lation function H{zi, Z2] Ax) together with its projections, is computed. It is demon- 
strated that the shape fluctuations follow the (extended) capillary-wave theory expres- 
sion S{q) = kT / {D + 'jq'^ + Kq'^) for sufficiently small wave vectors q. We find k negative, 
K < . At q = 2n there is a strong nearest-neighbor peak. Both these results confirm 
the earlier findings as established in simulations in d=3 and in continuous space, but 
now in an exactly soluble lattice model. 
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I.Introduction 

The shape fluctuations of interfaces are described by the capiUary-wave theory [1- 
10]. Its principal result is that the Fourier spectrum of the fluctuations follows for small 
wave- vectors the divergence. The divergence may be damped by an external fleld 
or by flnite size, as in the form 1/(D + 7g^), valid in the limit g ^ 0. Both D and the 
interfacial tension 7 are positive. The scattering experiments scan an interval of wave 
vectors up to the nearest-neighbor peak and beyond and many authors have searched 
for an extension of the simple capillary-wave form. An obvious generalization appeared 
to be 

S{q) = kT/{D + -fq^ + Kq^). (1.1) 

Sum rules require convergence of integrals over q and k > seems natural. Addi- 
tional input was provided by the fluctuations of other [d — 1)— dimensional objects in 
d— dimensional systems, such as membranes and bilayers in liquid solvents. According 
to both theory and experiment, the fluctuations of those autonomous objects are ruled 
not by 7 but by rigidity coefficient coupled to the curvature. It seemed only natural to 
"improve" the theory of fluctuations of ordinary non-autonomous interfaces by incor- 
porating a contribution of a hypothetical rigidity coefficient coupled to curvature. The 
latter contribution ought to be 

kT/nq'^ (1.2) 

with K necessarily positive. The positivity condition is necessary for the free energy of 
the deformed state to be larger than the free energy of the undeformed state. In this 
way one arrives at the expression (1.1) with the expectation of all three constants being 
positive. 

The first signs of inconsistency were signalled in a series of papers [8] where it 
was demonstrated that (for the long range attractive dispersion force, but not only) 
expression (1.1) is inadmissible and the last term must be replaced by nq^ f{q) where 
the function / contains nonalytical terms with logarithms. 

In a series of simulations by Molecular Dynamics of liquid-vapor surfaces, it was 
found independently[6] that all results pointed to a negative k. A negative value dis- 
proved the interpretation of k as a rigidity coefficient. The unexpected results were 
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recently confirmed by simulations of capillary waves in other liquid systems [9]. In ei- 
ther case[6,9] the intermolecular attractions between particles, were technically of short 
range because of the cutoflFs. 

It may be argued that these two independent simulation results [6, 9] are also con- 
firmed by scattering experiments; a negative k means that the denominator gets smaller 
than it would have been without any k, the scattering intensity gets bigger - and this was 
observed for liquid-vapor surfaces[10]. Such an interpretation, however, is not generally 
accepted and other, more complicated, explanations were sought [5, 10]. 

It would be interesting to confirm or negate the finding of a negative nhy a method 
other than simulation, for example in a system or model which is soluble exactly. 

II. The SOS model. 

In the Solid-on-Solid (SOS) model in two dimensions d=2 [11-17], a microscopic 
configuration of the interface is given by the collection of heights 

{hi},i = l,L; l<hi<M (2.1) 

as illustrated in Fig.l. The strip is M x L and below we consider the limit L ^ M. 
The other limit, M ^ has also been considered [14, 17]. The "vertical" variable is h 
or z at the horizontal position x = aoj, < j < L and we take ao = 1. The "height" is 
limited to < /i < M. The standard d = 2 SOS model[ll,17] assumes the energy of a 
given configuration as 

E ^ e^\hi+i - hiY" e>0 (2.2) 

with n = 1. A Gaussian model (n = 2) can also be considered. Introducing the M x M 
transfer matrix 

T(/i, /i') = gl'^-'^'l (2.3) 

with Q = exp[— e/A;T], we find that the partition function and correlation functions 
can be expressed in terms of ordinary scalar products of matrices T. If the transfer 
matrix is diagonalized - as it can be - all solutions are available, also as analytical 
expressions[ll] [14] [17] [15]. There is a minor and unimportant approximation involved 
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in the analytical expressions for the eigenvectors and eigenvalues of the inverse of T[15] 
(see Appendix). 

We have calculated the density-density correlation function H{zi, Z2; q) and other 
quantities such as the density profile p{z) and the probability functions: one-point p{h) 
and two-point P{hi, h2', Ax) and P{hi, h2', q). The following relations take place 

p{h) = p{h + 1) - p(/i); p{h) = Z{h, 1)2 (2.4) 

where Zi{h) = Z{h, 1) is the dominant eigenvector corresponding to the largest eigen- 
value Ai. Moreover 

H= J2 Yl P{hiM-p{hi)p{h2) (2.5) 

/ll>2l h2>Z2 

at each value of Ax = \x2 — xi\. Calculation of the Fourier Transform gives 

Hiz,,Z2;q)= E E/^(^)^^W^^(^2), (2.6) 

where 

CJ{h) = ZJ{h)Z^{h). (2.7) 

In 

fj{q) = (1 - r2)/(l + - 2r cos(g)), (2.8) 

r = Tj = Xj/Xi < 1 and we never need j = 1. Periodic boundary conditions in the 
direction were assumed, h{x = 1) = h{x = L + 1); next the limit L ^ 00 has been 
taken at constant M. It is useful for the discussion of the dependence to introduce 
projections which result in functions of one variable, q. One such function is H{q) which 
projects, at each q, the matrix H onto p{h). A similar projection was used in d=3[2,6]. 

^ = ^^^^(^1)^(^1' ^2)^/^2). (2.9) 

hi h2 

Since p{h) = Z\{K)^ , by using (2.5) we can change the order of summations to make a 
computationally convenient expression. 

4 



Another projection is the height-height correlation function < h{xi)h{x2) > whose 
Fourier transform is the well known structure factor (of the interface) 

S{q) =< hqh_q > . (2.10) 

We have 

< h{xi)h{x2) >=Y.Y1 hh'[P{h, h'; Ax) - p{h)p{h')] (2.11) 

h h' 

and, constructing the Fourier sums, 

h h' j>2 

where the largest eigenvalue is excluded from the sum over the eigenvalue index j. 
Alternatively the definition (2.6) of Cj can be invoked and the order of summation 
changed to give 

S{q)^J2fMdj{qf (2.13) 

i>2 

where 

d,(^)^^Z,(/.)/.Zi(/.). (2.14) 

h 

To summarize, we have so far the full density-density correation function H{zi, Z2,q), the 
corresponding derivative P{hi,h2,q) and the projections H{q) and S{q). Moreover, 
for a given value of q, the density-density correlation function ^isaM — IxM — 1 
positive definite matrix; it can be diagonalized and all its M — 1 eigenvalues are functions 
of q; the largest eigenvalue can be chosen and discussed as function of q. 

To make practical use of these expressions, the T— matrix is diagonalized as function 
of Q and this can be done analytically for small values of M, or one can take the route 
through the inverse of T, which is easier because the inverse is tridiagonal[ll]. The final 
expressions are given in the appendix[15]. 

III. Results for the wave-vector dependence of correlations. 

In order to reveal and illustrate the details of the capillary waves in the SOS model 
recalled and defined in the previous section, we numerically compute the analytical 
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expressions given in Section II and in the Appendix. In the Figures we show a selection 
of data, all for Q = 0.5. 

The numerical computation of H, of the projection H, and of the interfacial (string) 
structure factor S{q), reveals three regions of the variable q. These regions are clearly 
visible in Fig. 2. At the very lowest values of g a plateau (in logarithmic scale) practically 
equal to the limit for q = extends up to a crossover into the second and main 
region, that of the capillary waves. All the projections show a typical capillary-wave 
dependence of the form of (1.2), 1 / const. q^. This region extends up to ~ 1 — 5 where 
the amplitudes of the capillary waves have practically died out. Then in the third region 
{q^ > 5) the nearest neighbour peak appears near q ^2tt i.e. in Fig.2 near ~ 39 . 

The plateau and the finite value at (/ = 0+ arc due to the finiteness of the width 
M in the 2— direction perpendicular to the interface. The excursions of the interface 
from the average position in the middle at ~ M/2 are damped like waves in a channel, 
of width M and length L. This is captured by the smallest values of q. 

The main capillary-wave region dominates in the middle region and is the main 
physical effect. 

The wave- vector q is measured in units of inverse lattice constant ao, taken as unity 
here; hence the nearest neighbour peak expected at 27r/ao occurs at a; = ~ 477^. Its 
appearance was never considered in earlier work. 

This simple fact is, however, of importance because the presence of the peak requires 
a raise in the structure factor at g < 27r, and this gradual raise takes place already at 
values of q significantly lower than 27r. 

It may be found surprising that the raise begins so early. 

The plots of H for different sizes (widths) are shown in Fig. 3. The plot of S{q) in 
Fig.4 shows a curve common to all sizes in the central (capillary-wave) portion of the 
plot and the raise of S{q) preceding the n.n. peak. Each of the H{q) and S{q) for each 
size M were fitted to an empirical version of (1.1), i.e. to 

S{x) = \/{d + gx + K^x'^) x = q^ (3.1) 

where the constants cZ, or ^r, were estimated from the least-square fits. In all 
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cases but one the least-square fit produced a negative K4. The same formula was used 
for fitting the largest eigenvalue of H{zi, Z2', q) and invariably negative K4 was obtained. 
The raise in S(q) can be mistakenly attributed to other reasons, if the existence of the 
faraway nearest-neighbour peak is neglected. 

IV. Discussion 

We have found in this exactly soluble lattice model in d=2 dimensions the same 
feature discovered earlier in two[6,9] independent simulations: a raise in the structure 
factor over and above the capillary- wave contribution - which results in a negative value 
of the coefiicient k in (1.1), to which the empirical coeflBcient K4 is proportional . 
Therefore it cannot be interpreted as the rigidity of the interface. 

There exist systems with capillary waves and positive k, : as mentioned in the 
introduction these are bilayers (and membranes) in d=3 immersed in a liquid solvent. 
For each such autonomous object there exists a surface area which minimizes the free 
energy. No such minimum exists in non-autonomous interfaces nor in the Ising or SOS 
model. 

In general we must conclude that there exist several kinds of interfaces ruled by 
different laws. To strenghten this point one can invoke the existence of Langmuir mono- 
layers with their numerous surface phases. 

We may associate with the bilayers and membranes the term "autonomous" and 
with ordinary liquid- vapor and liquid-liquid interfaces the term "non-autonomous". 
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Appendix A. 

The T-matrix given by (2.3) produces on inversion a tridiagonal matrix r with 
diagonal elements given by 

d = (i + g2)/(i-g2) (Ai) 

and all off diagonal elements Tn,n+i,Tn,n-i given by 

o=-Q/(l-g2) (^.2) 

Strictly speaking not all diagonal elements are equal to d; the first and the last element 
are 

Til = TMM = do = 1/(1 - Q') (A.3) 

For large M, safely d can be substituted for do and the resulting matrix r' has the 
following eigenvalues [15] 

A,- = (1 + q2 _ 2gcos(6oi))/(l - Q') (A4) 

where bo = 7r/{N + 1) and j = 1,---,M. Hence j=l corresponds to the smallest 
eigenvalue, i.e. to the largest eigenvalue of the original matrix T, equal to 1/Ai. The 
eigenvectors are 

Zj{h) = {l/Af)cos{boju) (A5) 

for j odd , and 

Zj{h)^{l/X)sm{boju) (A6) 

for j even. These eigenvectors are orthonormal and complete. The dimensions of the 
matrices are M x M and M is best taken odd. The norm is Af^ = (M + l)/2 and 
i^ = h-l-{M + l)/2 so that = -(M + l)/2, 0, +(M+ l)/2. 

The eigenvectors are common to the diagonalized matrix and its inverse whereas 
the eigenvalues are the inverses. Because T is a symmetric positive definite matrix, it 
admits a representation 

Th,M = X^^,(/ii)(l/A,)^,(/i2) (A7) 
j 

and this representation for T was used to obtain expressions in Section II, remembering 
that (A. 4) refers to the inverse of T. 
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Figure Captions. 
Caption to Fig.l 

A portion of the SOS interface (which must be without overhangs or clusters). The 
lattice constant is 1, ordinates are x = j and abscissae are h{x) = hj < M. 

Caption to Fig.2 

An example of the projection H{q) plotted against x = in a, range wide enough 
to display all three regions; the plateau, the main capillary wave region H ~ 1/x; and 
the nearest neighbour peak at a; ~ 47r^. This particular curve for M — 31 was fitted to 
Eq.(3.1) with d = 0.2964, g = 122.40, K4 = -8.0. Q = 0.5 . 

Caption to Fig.3 

The projections H{q) plotted against x = (^ for the following widths M of the strip 
M= 120,001; 12,001; 1201; 601; 301; 121; 61; - from lowest values to the highest in that 
order. The line is j/ = l.e — 6/a; to which all B. are parallel, perfectly so in the center of 
the plot. All curves were fitted to eq.(l.l) and all fits produced negative ^4. Q = 0.5 . 

Caption to Fig. 4 

The interfacial structure factor Six) defined by (2.10) plotted against x = q^ for 
the widths M= 120001 (crosses), 12001 (plus signs), 1201 (stars), 121 (boxes), and the line 
fitted to (1.1). The nearest-neighbor peak and the faint increase in 5' for > 0.1 are to 
be noticed. This faint increase forces a negative k. Q = 0.1. The slope is practically 
equal to the stiffness coefficient [17], which explicitely is (1 — Q'^)/{2Q). 
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